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Abstract 

y 7 " 1 Geometrical and topological properties of the quasi-plane Szekeres model and of 

^Jy the plane symmetric dust model are discussed. Some related comments on the quasi- 

hyperbolical model are made. These properties include: (1) The pattern of expansion 
in the plane symmetric case, and the Newtonian model that imitates it; (2) The 
possibility of toroidal topology of the t = const sections in the plane symmetric model; 
(3) The absence of apparent horizons in the quasi-plane and quasi-hyperbolic models 
(they are globally trapped); (4) Description of the toroidal topology in the Szekeres 
coordinates; (5) Consequences of toroidal topology in the nonsymmetric quasi-plane 
. model; (6) Avoidance of shell crossings in the toroidal model; (7) Interpretation of 

the mass function in the quasi-plane model, with the toroidal and with the infinite 
^ ' space. 
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X 1 Motivation 



The quasi-spherical Szekeres model [T] - [T7] is rather well-understood by now. In spite of 
its nontrivial geometry, its basic defining features are not too difficult to grasp intuitively. 
In a simple-minded way one may say that it is obtained when the spherical symmetry orbits 
in the Lemaitre - Tolman model [HI [T7j are made nonconcentric to destroy the symmetry, 
but the energy- momentum tensor is still that of dust. Recently, that model even found 
application to solving problems directly related to observational cosmology [US, E]- in 
contrast to this, the first serious attempt to interpret the quasi-plane and quasi-hyperbolic 
models [18] revealed that even the corresponding plane- and hyperbolically symmetric 
models are not really understood and require more investigation. Some properties of those 
models were established in Ref. [18], the present paper is a continuation of that research. 

The aim of the present paper is to clarify some of the basic geometrical features of the 
quasi-plane Szekeres model, and of the plane symmetric dust model. The following topics 
are investigated here: (1) The pattern of expansion in the plane symmetric model, and 
the Newtonian model that imitates it; (2) The possibility of toroidal topology of the t = 
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const sections in the plane symmetric model; (3) The absence of apparent horizons in the 
quasi-plane and quasi-hyperbolic models (they are globally trapped); (4) Description of 
the toroidal topology in the Szekeres coordinates; (5) Consequences of toroidal topology 
in the nonsymmetric quasi-plane model - it turns out that in the nonsymmetric case the 
orbits of quasi-symmetry must be projective planes rather than tori; (6) Avoidance of 
shell crossings - in the toroidal model they can be avoided, unlike in the model with 
infinite orbits; (7) Interpretation of the mass function in the quasi-plane model. For the 
most part, the paper is devoted to showing that the space of constant time in the plane 
symmetric dust and the quasi-plane Szekeres models can be interpreted as a family of flat 
tori, with the ones of smaller diameter enclosed inside those of larger diameter. Such a 
topology explains several properties of the models, among them the pattern of accelerated 
expansion and the finiteness of the mass function, and also it allows to construct models 
that are free of shell crossings. In the end it turns out that these models are of lesser 
use in astrophysical cosmology than the quasi-spherical ones. Because of being globally 
trapped, they cannot be used for modeling dynamical black holes. Because they expand 
by the same law as the positive-energy Lemaitre - Tolman model, they cannot model the 
formation of structures that collapse to very dense states. They might be applicable for 
the description of formation of moderate condensations, like galaxy clusters, and of voids. 

Mena, Natario and Tod have recently also considered the quasi-plane and quasi-hyper- 
bolic Szekeres models with toroidal and higher- genus topologies [Tj5]. They considered the 
matching of those solutions, with nonzero cosmological constant (corresponding to A > 
in the notation adopted here), to the plane- and hyperbolically symmetric counterparts 
of the Schwarzschild solution, also allowed to have nontrivial topologies of the symmetry 
orbits. However, there is no overlap between their results and those of the present paper, 
as they mainly considered the global geometry of the resulting black hole, while here the 
emphasis is put on local geometry of the topologically nontrivial Szekeres spacetime. 

In the next section, basic facts about the Szekeres solutions are recalled. 

2 Introducing the Szekeres solutions 

The metric of the Szekeres solutions is 

ds 2 = dt 2 - e 2a dz 2 - e 2f3 (dx 2 + dy 2 ) , (2.1) 

where a and (3 are functions of (t, x, y, z) to be determined from the Einstein equations 
with a dust source. The coordinates of (12.1 j) are comoving so the velocity field of the dust 
is u M = 5^0 and u M = 0. 

There are two families of Szekeres solutions, depending on whether /3, z = or f3, z ^ 0. 
The first family is a simultaneous generalisation of the Friedmann and Kantowski - Sachs 
[2"U] models. So far it has found no useful application in astrophysical cosmology, and we 
shall not discuss it here (see Ref. [IS]); we shall deal only with the second family. After 
the Einstein equations are solved, the metric functions in (12. ip become 

e a = h(z)$(t,z)P, t =h(z)($, t +$v, t ), 
e' u = A{z) (x 2 + y 2 ) + 2B 1 (z)x + 2B 2 (z)y + C(z), (2.2) 
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where the function $(£, z) is a solution of the equation 



2 , , , 2M(z) 1 , n , N 

while /i(z)? M(z), A(z), Bi(z), B 2 (z) and C(z) are arbitrary functions obeying 

g(z) d ^ AC - £? x 2 - B 2 2 = I [l/tfO*) + *(*)] . (2.4) 

The mass- density is 

(2Me 3 ^ 



In the present paper we will mostly consider the case A = 0. 

This solution has in general no symmetry, and acquires a 3-dimensional symmetry 
group with 2-dimensional orbits when A, B\, B 2 and C are all constant (that is, when 
u, z = 0). The sign of g(z) determines the geometry of the surfaces (t = const, z = const), 
and the symmetry of the limiting solution. The geometry is spherical, plane or hyperbolic 
when g > 0, g = or g < 0, respectively. With A, B 1 , B 2 and C being functions of z, the 
surfaces z = const within a single space t = const may have different geometries (i.e. they 
can be spheres in one part of the space and the surfaces of constant negative curvature 
elsewhere, the curvature being zero at the boundary). The sign of k(z) determines the 
type of evolution; with k > = A the model expands away from an initial singularity and 
then recollapses to a final singularity, with k < = A the model is either ever-expanding 
or ever-collapsing, depending on the initial conditions; k = = A is the intermediate case 
corresponding to the 'flat' Friedmann model. 

The Robertson- Walker limit follows when z = r, <&(t,z) = rR(t), k = k$r 2 where 
ko = const and B\ = B 2 = 0, C = 4 A = 1. This definition of the R-W limit includes 
the definition of the limiting radial coordinate (the Szekeres model is covariant with the 
transformations z = f(z'), where f(z') is an arbitrary function). 

The Szekeres models are subdivided according to the sign of g(z) into the quasi-spherical 
ones (with g > 0), quasi-plane (g = 0) and quasi- hyperbolic (g < 0)0 Despite suggestions 
to the contrary made in the literature, the geometry of the latter two classes has, until 
very recently, not been investigated at all and is not really understood; work on their 
interpretation has only been begun by Helalby and Krasinski [18]. The sign of g(z) is 
independent of the sign of k(z), but limitations are imposed on k(z) by the signature of 

the spacetime: for the signature to be the physical (H ), the function h 2 must be 

non-negative (possibly zero at isolated points, but not in open subsets), which, via ( 12. 4p 
means that g(z) — k(z) > everywhere. Thus, with g > (in the quasi-spherical case) 
all three possibilities for k are allowed; with g = only the two k < evolutions are 
admissible, and with g < 0, only the k < evolution is allowed. 

Only the quasi-spherical model is rather well investigated, and found useful application 
in astrophysical cosmology. We recall now its basic properties. 



^^We stress once again that the same Szekeres model may be quasi-spherical in one part of the spacetime, 
and quasi-hyperbolic elsewhere, with the boundary between these two regions being quasi-plane, see an 
explicit simple example in Ref. [18 . In most of the literature published so far, these models have been 
considered separately, but this was only for purposes of systematic research. 
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It may be imagined as such a generalisation of the Lemaiitre - Tolman (L-T) model in 
which the spheres of constant mass were made non-concentric. The functions A[z), Bi(z) 
and B 2 (z) determine how the center of a sphere changes its position in a space t = const 
when the radius of the sphere is increased or decreased (see a more detailed discussion 
of this in Ref. [15J). Still, this is a rather simple geometry because all the arbitrary 
functions depend on just one variable, z. They give us some limited possibility to model 
the real structures in the Universe (see very elegant examples in Refs. [131 ED; but a 
fully satisfactory model should involve arbitrary functions of all three spatial variables, to 
allow modelling of arbitrary structures. Such models are still nonexistent, so the Szekeres 
models are so far the best devices that exist. 

Often, it is more practical to reparametrise the arbitrary functions in the Szekeres 
metric as follows (this parametrisation was invented by C. Hellaby [21]). Even if A = 
initially, a transformation of the (x, y)-coordinates can restore A ^ 0, so we may assume 
A 7^ with no loss of generality (see Ref. [16J). Then let g ^ 0. Writing 



(A,B 1 ,B 2 ) 



\9\ 



S 



-P,-Q), 



def , I I 

£ = 9/\9U 



k = -4\g\ x 2E, 
we can represent the metric ( 12. ID as 



p def S 

t ~2 



ds 2 



(R, z —RS, z /£) 



M = 8\g\ 3/2 M and $ = Ry/\4g\, (2.6) 



x 



P 



s 



+ 



y-Q 
s 



+ e 



dt 



R 



e + 2E(z) 



dz 2 - [dx 2 + dy 2 ) . 



[2.7) 



When g = 0, the transition from fl2H) to (J32D is A = 1/(2S), B x = -P/(2S), B 2 = 
-Q/(2S), k = -2E, M = M and $ = R. Then (ETJ) applies with e = 0, and the resulting 
model is quasi-plane. 

For further reference, the evolution equation (|2.3|) . in the variables of (12.71) becomes 



, , N 2M(z) 1 , 9 
R, t 2 = 2E(z) + — ±1 + -AR 2 ; 



(2- 



From now on, we will use this representation; the tilde over M will be dropped, but it 
must be remembered that the M in ( 12.81) is not the same as the one in ( 12. 3ft . 

The representation (12.71) makes the calculations simpler because the arbitrary functions 
in it are independent (the condition (12.41) has been incorporated in this form). However, 
it creates the illusion that the cases e = +1,0,-1 characterise the whole spacetime, and 
thus obscures the fact that all three cases can be parts of the same spacetime. 

Rotation and acceleration of the dust source are zero, the expansion is 



e 



^R,t R,tz —R,t R,z I R 
R R, z —R£,z/£ 



(2.9) 



and the shear tensor is 



a a p = -S diag (0,2,-1,-1), where 
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$, t2 -$, t $, 2 /$ ^ R, tz -R, t R, z /R 

R, z -R£, z /S ' 1 ' j 

Definitions of the Szekeres solutions by invariant properties can be found in Ref. [TO] . 

When A ^ 0, the solutions of (12. 8p involve elliptic functions. A general formal integral 
of (12.81) was presented by Barrow and Stein-Schabes [22]. Any solution of (12.81) will contain 
one more arbitrary function of z that will be denoted ts(^), and will enter the solution in 
the combination (t — tsiz)). The instant t = te(z) defines the initial moment of evolution; 
when A = it is necessarily a singularity corresponding to $ = 0, and it goes over into 
the Big Bang singularity in the Friedmann limit. When t^, z ^ (that is, in general) the 
instant of singularity is position-dependent, just as it was in the L-T model. 

Just as in the L-T model, another singularity may occur where (e^) , 2 = (if this 
equation has solutions). This is a shell crossing, but it is qualitatively different from that 
in the L-T model. As can be seen from (12.21) . in the quasi-spherical case, when a shell 
crossing exists, its intersection with a t = const space will be a circle, or, in exceptional 
cases, a single point, not a sphere. In the quasi-spherical models shell crossings can be 
avoided altogether if the arbitrary functions are chosen appropriately, see the complete list 
and derivation in Ref. |15j . In the quasi-hyperbolic models, shell crossings can be avoided 
in one sheet of each hyperboloid, but are unavoidable in the other, see Ref. [18]. In in 
the quasi-plane model, if the flat surfaces existing in it are interpreted as infinite planes, 
shell crossings are unavoidable [18]. However, if these surfaces are interpreted as being 
compact, shell crossings can be avoided - see Sec. [9] in the present paper. 

Equation (12.81) is identical with the Friedmann equation, but, just like in the L-T limit, 
with k and M depending on z, each surface z = const evolves independently of the others. 

The models defined by (12.11) - (12. 5p contain 8 functions of z, but only 5 of them corre- 
spond to independent physical degrees of freedom. One of the 8 functions is determined 
by (12. 4p . g(z) was made constant by the reparametrisation (12.61) . and one can be specified 
by a choice of z, for example by defining z' = M, or M = z' 3 x {a constant}. 

A quasi-spherical Szekeres region can be matched to the Schwarzschild solution across 
a, z = const hypersurface [6]. The other two Szekeres regions can be matched to the plane- 
and hyperbolically symmetric counterparts of the Schwarzschild solution (see Ref. [23] for 
the solutions and [18] for the matching). 

In the following, we will represent the Szekeres solutions with /3, z ^ in the parametri- 
sation introduced in (12.71) . The formula for density in these variables is 

_ 2 (M, z -3MS, Z IS) 
Kp ~ RHR, z -R£, z /£)> (2 - 11} 

where, let it be recalled, the M above is the M of (12. 6ft . 



3 The plane symmetric models 

The plane symmetric dust models (first found by Ellis [21]) result from (12.71) when e = 
and (P, Q, S) are independent of z. The constant S can then be scaled to 1 by appropriate 
redefinitions of R, E and M. Then, with constant P and Q, the coordinate transformation 

x = P+-^— 2 , y = Q + ^L- (3.1) 
p z + q p + q 
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changes the metric to 



ds 2 = dt 2 - ^yd^ 2 - R 2 (dp 2 + dq 2 ) , (3.2) 



while the energy-density simplifies to 

— p = Wr7 z - (3 ' 3) 

These models are called plane symmetric because their symmetries are the same as those 
of the Euclidean plane; in the coordinates of ( 13. 2ft they are: 

p' = p + A u (3.4a) 
q = q + A 2 , (3.4b) 
ip',q') = (pcosa + gsina, — psina + gcosa), (3.4c) 

where Ai, A 2 and a are arbitrary constants - the group parameters. 

Note that eqs. (12.81) and (I3.3P are identical to their counterparts in the spherically 
symmetric models. In particular, the function M(z) enters in the same way as the active 
gravitational mass did in spherical models. However, if we wish to interpret M(z) as a 
mass contained in a volume, we encounter a problem - see below. 

Examples of plane symmetric spaces are the Euclidean plane and the Euclidean space 
i?3 with the metric ds3 2 = dx 2 + dy 2 + dz 2 . However, the space of constant t in ( 13.21) can 
never become flat; its curvature tensor is [18J: 

3r>zp _3ra 3 r>pq _ (nr\ 

Jl zp — 11 zq — j^j^ j Xh pq — ^ 2 . \°-°) 

Nevertheless, the surfaces P 2 of constant t and z in ( I3.2p are flat. Thus, there is some 
mystery in the geometry of the spacetimes (13.21) . One component of the mystery is this: In 
the quasi-spherical case, and in the associated spherically symmetric model, the surfaces 
of constant t and z were spheres, and M(z) was a mass inside a sphere of coordinate radius 
z. In the plane symmetric case, if the P 2 surfaces are infinite planes, they do not enclose 
any finite volume, so where does the mass M(z) reside? 

With M being a constant, the metric (13. 2p becomes vacuum - the plane symmetric 
analogue of the Schwarzschild spacetime. 

In the quasi-spherical Szekeres, and in spherically symmetric, solutions, analogies exist 
between the relativistic and the Newtonian models. We will thus compare the plane 
symmetric model with its possible Newtonian counterparts. For this purpose, let us note 
the pattern of expansion in (13. 2p and (12.81) with A = 0. When R, t ^ 0, eq. (I2.8P implies 

R, tt =-M/R 2 . (3.6) 

Note that R,u = implies M — 0, which is a vacuum (in fact, Minkowski in strange 
coordinates) metric. Now take a pair of dust particles, located at (t, zi,p ,q ) and at 
(t, z 2 ,p ,q ), and consider the affine distance between them: 

Mt)= r^=,^ = (3.7) 



: , v2E dt 2 J zi ^2E 
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Thus, the two particles will be receding from each other (or approaching each other if 
collapse is considered) with acceleration that can never be zeroo 

Take another pair of dust particles, located at (t, z ,pi,q ) and at (t, £o>P2, Qo)- The 
distance between them, measured within the symmetry orbit, is 

Rdp = R(p 2 - Pi) — ^ = R, tt (P2 - Pi) , (3.8) 

i.e. the acceleration of the expansion can never vanish in this direction, either, unless 
M — 0. The same result will follow for any direction in the (p, q) surface. Thus, the 
expansion or collapse in this model proceeds with acceleration in every spatial direction. 
We will compare this result with the Newtonian situation. 



4 A Newtonian analogue of the plane symmetric dust 
spacetime 

At first sight, it seems obvious that the Newtonian model analogous to the plane symmetric 
dust model in relativity should be dust whose density is constant on parallel (x, y)-planes, 
and depends only on z. Let us follow this idea. 

If the potential is plane symmetric, then, in the adapted coordinates, it depends only 
on z. Thus, the Poisson equation simplifies to 

d 2 V 

-^ = 4ttG P (z). (4.1) 

The general solution of this is 

V = AuG I dz I dzp{z) + Az + B, (4.2) 

J z J z 

where A and B are integration constants; Zq is a reference value of z at which we can 
specify an initial condition. If we wish to have V = const (i.e. zero force) when p = 0, we 
must take A = 0, and then V(zq) = The equations of motion in this potential are 

dv* dV . * 

where v 1 are components of the velocity field of matter, so 

dvX dvV n dvZ dV a r~\ f Z J ' / '\ 

— - — = — - — = 0, — = — — = -AttG / dzp(z'). (4.4) 
dt dt dt dz J Z0 Iy 1 y ' 

This, however, gives a pattern of expansion different from that in the relativistic plane sym- 
metric model. In the above, expansion with acceleration proceeds only in the ^-direction, 
while in the directions orthogonal to z there is no acceleration, or, in a special case, not 



2 The acceleration would be zero if R,ttz — 0, which leads to a contradiction in the Einstein equations. 
3 An A ^ would be qualitatively similar to the cosmological constant in relativity. 
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even any expansion. Consequently, no obvious Newtonian analogue exists for the relativis- 
tic plane symmetric modeljj 

Equation (14. 4p shares one property with the relativistic evolution equation (12.81) . If p 
is bounded in the range of integration, then the force that drives the motion of the fluid is 
finite, giving the illusion that the potential is generated by some finite mass. However, if we 
wanted to calculate V(z) by summing up contributions to it from all the volume elements 
of the fluid, like is done in calculating the gravitational potentials of finite portions of 
matter, then the result would be an infinite value of V, in consequence of the source 
having infinite extent in the (x, y)-plane. Thus, if we want to interpret the r.h.s. in (14. 3p 
as being generated by a mass, then the mass that drives the evolution is not the total mass 
in the source, but the mass of a finite portion of the source. 

We now provide a solution of the Poisson equation that qualitatively mimics the pattern 
of expansion of the plane symmetric relativistic model. The equipotential surfaces of the 
potential in question will be locally plane symmetric, but their symmetries will not be 
symmetries of the whole space. 

Consider two families of cones given by the equations (see Fig. [I]) 



// : z nr. r =z + r/a, r = f \/ x 2 + y 2 , (4.5) 



where a is a constant and (x, y, z) are Cartesian coordinates. The cones of constant u 
are orthogonal to the cones of constant v , and the two families are co-axial. We choose 
u and v as two coordinates in space; the third coordinate will be the angle ip around 
the axis of symmetry. We begin with the Euclidean metric in the cylindrical coordinates, 
ds 2 = dr 2 + r 2 d<p 2 + dz 2 , and transform this to the (u, v, <p) coordinates by 



a(v — u) u + a 2 v 



1 + a 2 ' " "~ 1 + a 2 



(4.6) 



The transformed metric is 



= du 2 + a 2 dv 2 a 2 {v-u) 2 dp 2 

1 + a 2 (1 + a 2 ) 2 ' 1 ' ' 



The Laplace operator, which in the cylindrical coordinates is 



. I d ( dV\ 1 d 2 V d 2 V 

AV = r H 1 , (4. 

r dr \ dr J r 2 dp 2 dz 2 



in the (u, v , (/^-coordinates becomes 
AV = (1 + a 2 ) 



1 dV d 2 V 1 dV 1 d 2 V 1 + a 2 d 2 V 

v — u du du 2 a 2 {v — u) dv a 2 dv 2 a 2 (v — u) 2 dp 2 



Thus, if V depends only on u, then the Poisson equation says: 



(4.9) 



(4,0) 



V 



4 Incidcntally, there will be no Newtonian analogue for the hyperbolic model, since the orbits of hy- 
perbolic symmetry cannot be embedded in a Euclidean space at all. They can be embedded in a flat 
3-dimensional space, but the space then must have the signature ( — h +). 
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Fi gure 1: The cones where u is constant are orthogonal to the cones where v is constant. The coordinates 
in space are u, v and the angle around the axis of symmetry. The figure shows an axial cross-section 
through the setup. A gravitational potential which depends only on u in these coordinates gives an 
expansion pattern that is qualitatively similar to the one in a plane symmetric dust spacetime. 



This density has a singularity at u — v, i.e. at r = in the cylindrical coordinates. 

The gradient of V(u) has nonzero components in all directions, and so will create 
expansion decelerated in all directions. The expansion will be isotropic with respect to the 
u = v axis, and the anisotropy between the (x,y) and the z-directions is controlled by a. 

This potential was introduced here for illustrative purposes, and it is doubtful that it 
could correspond to any physical configuration. In order to make it credible, one should 
solve the continuity equation and the Euler equations of motion in it. We do not quote 
here the appropriate calculations because they lead to an intransparent tangle of differential 
equations. For dust, that set of equations is overdetermined, so probably has no solutions. 

5 Plane symmetric 3-spaces interpreted as tori 

Although known for a long time (see Ref. [24J), the plane symmetric model has not been 
investigated for its geometrical and physical properties. 

Since a flat spatial geometry is not possible in it (see eq. (13. 5p ). we now consider other 
possible 3-geometries with planar symmetry. The next simplest is a space of constant 
curvature. From ( 13. 51) . the space of constant t = to will have constant curvature when 

2E = ±C 2 R 2 , (5.1) 

where C is a constant. The curvature is positive when E < and negative when E > 0. 
Since the signature of spacetime requires E > 0, we follow only the + case. Choosing 



9 



R(t , z) = R as the spatial coordinate in this space, we get: 



dsS = S 2 



2 



dR 
C 2 R 2 



+ R 2 (dp 2 + dq 2 ) 



(5.2) 



Note that only one hypersurface can have the property (15. ip (since E is independent of t 
while R depends on t). Thus, the 3-geometry of a space of constant t can evolve away from 
or toward (15.21) . or through (15. 21) . but cannot preserve this geometry over a finite time. 

The surface of constant q in ( 15. 2ft has the metric ds2 2 = di? 2 / (C 2 R 2 ) + R 2 dp 2 . To 
visualise it, we embed it now in a 3-dimensional Euclidean space with the metric 

ds 3 2 = dZ 2 + di? 2 + R 2 dp 2 . (5.3) 

Our ds 2 2 is the metric of the surface Z = Z (R), where Z 0tR 2 + 1 = l/(Ci?) 2 , thus 



yi ^cm i 



CR 



In + Vl^C^R 2 

1 + VI - C 2 i? 2 



(5.4) 



This embedding is possible only in the range R < 1/C. The R > 1 jC part of the surface 
can be embedded in a flat space of the signature ( — h +)Jj The metric is then 

ds 3 2 = -dZ x 2 + di? 2 + R 2 dp 2 , (5.5) 

and the embedding equation is 



Z 1 = ± 



/\/CJ 2 R 2 — 1 1 r / \ 

V CR di? = ±- [VC 2 i? 2 - 1 - 2 arctan [CR + VC 2 i? 2 - lj + tt/2 

(5.6) 

(the constant of integration was chosen so that Z Q (1/C) = Z\{l/C)). The functions Z(R) 
and Zi(R) are shown in Fig. [21 Note that in both embeddings, (15.31) and (15.51) . p appears 
as the polar angle in the plane (R,p). If p is to be interpreted as actually being a polar 
angle, with the period 2n, then all points with the coordinates (t,z,p + 2irn,q), where n 
is any integer, should be identical with the point of coordinates (t, z,p, q). Since p — > p+ 
constant are symmetry transformations of the spacetime (3.1), there is no problem with 
such an identification. Thus we should imagine the (R,p) surface as being created by 
rotating the curve from Fig. [2] around the Z axis. 

However, the same picture would be obtained for an (R,q) surface in (15.21) . given by 
p = const. We would find that in that surface, q is the angular coordinate of the polar 
coordinates (R,q), and points of coordinates (t,z,p,q + 2nm) can be identified with the 
point of coordinates (t, z,p, q). We are thus led to conclude that (p, q) are both angular co- 
ordinates with the period 2n, and that the points of coordinates (p, q) have to be identified 
with the points of coordinates (p + 2nir, q + 2rmr), where n and m are arbitrary integers. 
The tentative conclusion is that the (p, g)-surface is a flat torus. 

The conclusion is tentative in the sense that, while we identify the set p = Po with the 
set p = po + 2tt, we are still free to carry out the symmetry transformations within the set 
p = pq. Thus, the identification can possibly be done with a twist, that will turn a square 



5 A similar phenomenon is known from the maximally extended Reissner-Nordstrom spacetime, when 
the region inside the interior horizon is depicted, see Ref. 16J. 
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Figure 2: Left panel: The function Z (R) (from to 1, solid line) and the function Z\{R) (from 1 to 
3, dotted line), given by eqs. (|5.4p and (|5.6p , respectively. In the graph we chose the — sign for Z (R) and 
the + sign for Zi(R). The graph is the cross-section of the (i?,p)-surface in the spaces (|5.3[) (left part of 
the curve) and (|5.5p (right part of the curve). Right panel: The (R,p) (or (R,qj) surface obtained by 
rotating the graph from the left panel around the R — axis. The lower end of the funnel is where the 
embedding in the Euclidean space breaks down, i.e. where the solid line meets the dotted line in the left 
panel. Upwards, the funnel goes infinitely far and becomes infinitely thin. 



into a Mobius strip, or with a two-way twist, that will turn it into a projective plane. We 
will see later that the latter case is in fact forced upon us in the full nonsymmetric case. 
We will use the term "toroidal topology" that will be meant to include an ordinary torus, 
and also the identifications with twists. 

The conclusions drawn from an embedding can be misleading. As an example, consider 
the hyperbolically symmetric counterpart of (13.21) . which can be written in the form: 

ds 2 = dt 2 - ^jL_ ds a - R 2 (dtf 2 + sinh 2 Mtp 2 ) . (5.7) 

The surface of constant t and constant cp has the metric ds2 2 = R, z 2 dz 2 / (2E(z) — 1) — 
R 2 di9 2 , and embedding it in a Euclidean space we would conclude that $ is the polar 
coordinate with the period 2tt. However, in this case $ — > i?+ constant are not symmetry 
transformations of the spacetime (or of a space of constant t), and so identifications of 
points with different values of •& are not permitted. Thus, the embedding in this case is not 
a one-to-one mapping. Consequently, the toroidal interpretation of the plane symmetric 
case must be treated as one possibility, and not as a definitive conclusion. 

Note from (15. 2p that the length of any segment of a curve given by p = const and q = 
const that goes into the point R = is infinite. The length of such a "radial" line between 
the values Ri and R 2 is 



f R2 dR 




Hf) 




J to CR 









Thus, Fig. [2] correctly suggests that a surface of constant p or constant q in the space (15. 2ft 
has the shape of an infinite funnel, and the point with coordinate R = is not accessible 
(does not in fact belong to this surface). This conclusion is consistent with the observation 
made in Ref. [T$J that in the planar Szekeres metric "there is no real origin, but R, M 
and E can asymptotically approach zero'fl 

6 In this quote, notation has been adapted to the presently used one. 
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Figure 3: A sketch of the full 3-space of the plane symmetric toroidal model (a faithful picture cannot be 
drawn because the 2-dimensional flat torus cannot be embedded in a Euclidean space, and the 3-space of a 
planar Szekeres model cannot be made flat). Each square section of the funnel represents a 2-dimensional 
flat torus, so its front edge coincides with the back edge, and the left edge coincides with the right edge. 
Each torus is an orbit of the symmetry group of the model. Each smaller torus is contained within all 
the larger ones, but the inclusion relation cannot be depicted in the Euclidean space of the picture. Also, 
in the curved 3-space, the 2-tori surround the asymptotic origin, which is the tip of the funnel, situated 
infinitely high above the plane drawn in the picture. 



It can be seen from (13. 2p that the (p, q) surface should have a toroidal topology with 
any form of E. The 3-metric of a t = to space is: 

d, 3 2 = + Ro 2 (dp 2 + dg 2 ) , (5.8) 

where -Ro(^) == R{to, z )- We can now embed a surface of constant p or a surface of constant 
q in a 3-dimensional flat space by the same method that we used for (15.21) . only the equation 
of embedding will not be explicit, it will be 

± Zj + 1 = 1/{2E) > 0, (5.9) 

the upper sign being for embedding in the Euclidean space, the lower sign for the embedding 
in the pseudoeuclidean space. In each case the coordinates p and q turn out to be the 
azimuthal coordinates. As argued in Ref. [18], if a nonsingular origin (where R = E = 0) 
is to exist, then it will be infinitely far from every point of the t = const space. This 
implies the infinite funnel geometry of Fig. [2J A sketch of such a torus is shown in Fig. [3J 
The toroidal geometry and topology of the (p, q) surfaces neatly explains the pattern 
of expansion. The circumference of the torus along the p- or g-direction is 2nR in the co- 
ordinates of (15.81) . Thus, as R increases with time, the circumference of the torus increases 
in proportion to R, which causes transversal expansion. 
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The toroidal topology also solves the problem of where the mass generating the gravita- 
tional field resides. As observed in Ref . [18] , the regularity conditions at an origin, derived 
for the Szekeres models in Ref. [15], are independent of e (it is not guaranteed that they 
can always be fulfilled; with e < they cannot). Thus, E/M 2 / 3 must tend to a nonzero 
constant as z — > z c , where z c is the origin. Knowing this, let us calculate the amount of rest 
mass in an arbitrary volume V, from (I3.3P and (13.21) . That amount is M. == j v p^/g^d^x, 
where g 3 is the determinant of the 3- metric of a t = const subspace of (13.21) . Thus 

c 2 r m c 2 r 1 

M = — — / -^Ldpdqdz = — — / ^=dpdqdM. (5.10) 



47rG J v ^f2E 4:wG J v y/2E 

With E oc M 2//3 in the vicinity of M = 0, the integral with respect to M is finite. With a 
toroidal topology, the ranges of p and q are finite, so the integrals over p and q also give 
final values. Thus, the total amount of mass in each space t = const is finite. 

The relation M. 1Z = M, z / V2E that follows from ( 15.101) is similar to the relation M., T = 
M, r I sj\ + 2E that held in the spherically symmetric and quasi-spherical models. By 
analogy with that case, we may thus say that in the plane symmetric spacetime the factor 
l/y2E measures the relativistic mass defect/excess, i.e. the discrepancy between the 
active gravitational mass M and the sum of rest masses A4. 

6 No apparent horizons in the quasi-plane and quasi- 
hyperbolic models 

An apparent horizon is the envelope of the region of trapped surfaces. A (past or fu- 
ture) trapped surface is such a surface on which both the ingoing and outgoing (past- 
or future-directed respectively) families of null geodesies converge. A future AH always 
forms in spherically symmetric or quasi-spherical-Szekeres collapse before the Big Crunch 
singularity is achieved, a past AH always exists after the Big Bang singularity. 

It turns out that the AH-s do not exist in the quasi-plane and quasi-hyperbolic Szekeres 
models, and, consequently, neither do they exist in the plane- and hyperbolically symmetric 
dust models. Actually, a stronger result holds: these spacetimes remain trapped for all 
the time. This follows by the method invented by Szekeres [2], which applies here almost 
unchanged - only the final conclusion is radically different in consequence of the different 
sign of e. To avoid getting into complicated details, we begin by using the general form 
(12.11) of the metric. Suppose a trapped surface exists, and call it S. 

We assume £ to be one of the orbits of the quasi-symmetry, i.e. to have its equation 
of the form {t = constant, z = constant}. It will be explained later (see after (16. 9p ) 
why it is sufficient to consider such surfaces to prove the conclusion. The traditional 
definition of a trapped surface requires that it be compact. With the toroidal topology in 
the planar model, our £ will be compact indeed. With the infinite topology, and in the 
quasi-hyperbolic model, E will be infinite. In view of the final result of our consideration, 
this fact will turn out to be unimportant. We choose these infinite surfaces because of 
their simple geometry. 

Consider any family of null geodesies intersecting S orthogonally, and let the tangent 
vector field of those geodesies be . Let (t,z,x,y) = (x°, x 1 , x 2 , x 3 ). We have then 

fc M F = 0, fcT; F =0 everywhere (6.1) 
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because k^ is tangent to null geodesies, and 

e = e = 0} (k ) 2 - e 2a (k 1 ) 2 = onE (6.2) 

because fc M is assumed orthogonal to E, so at the points of E it must be spanned on the 
vector fields normal to S, which are (1,0,0,0) and (0,e _a ,0,0). The affine parameter 
along each null geodesic may be chosen so that 

k° = e Q , k 1 = e = ±1 on E, (6.3) 

where we will call the geodesies with e = — 1 "ingoing" , and those with e = +1 "outgoing" 
A surface E is trapped when the expansion k^;^ calculated on E is negative for both 
families. We have on E, using (16.21) : 

fc% = k°, t +k\ z +e a (a, t +2{3, t ) + e (a, z +2/5,, ) . (6.4) 

In order to simplify this, we now differentiate the first of (16.11) by t, and write out the 
second of (16. ip for \i = 1, in both cases taking the result on E, i.e. making use of the 
simplifications given in (16. 2\i : 

k°, t -ee a k\ t -e a a, t =0, 

e a k\ t +e (k\ z +2e a a, t )+a„= 0. (6.5) 
Eliminating k l )t from (16. 5p . and using the result to substitute for k°, t +fc 1 ,z in (16.41) we get 

k^=2(e a /3, t +e(3, z ). (6.6) 
Using now the expressions for e a and e 13 in the notation of (12.71) . i.e. 

gQ = R, z -R£, z /£ ^ = R 

Ve + 2E(z) ' £ 1 ' ; 

we get in (16.61) 

The first factor changes sign only at shell crossings (see Ref. [15]), so we take it to be 
positive. Consider collapse, R, t < 0. For the ingoing family, e = —1, we have & M ; M < 0, 
without further conditions. For the outgoing family, e = +1, k^]^ will be negative when 
R,t j\Je + 2E < —1, which, with negative R, t , means that R, t 2 > e + 2E. Using (12.81) with 
A = for i?,t 2 , we then obtain 

2M/R > e. (6.9) 

With e = and e = —1, this is always fulfilled^] with the only exception of the 'asymptotic 
origin' in the planar model, where M/R = e = 0. 



7 When the surface of constant t and z is infinite, it cannot be closed, therefore the labelling "ingoing" 
and "outgoing" is only conventional. 

8 Note that M must be positive, or else (|3.6[) would imply that collapse is retarded and expansion 
accelerated. This would be gravitational repulsion. 
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A surface given by {t = constant, z = constant} passes through every point of the 
spacetime. Since each such surface has now been shown to be trapped at all of its points, 
this means that all points of the whole spacetime are trapped. 

Thus, the quasi-hyperbolic and quasi-plane model, along with their hyperbolically- and 
plane-symmetric limits, are globally future-trapped (when collapsing), and no apparent 
horizon exists for them. It follows now easily that the corresponding expanding models 
are globally past-trapped. 

This is consistent with the fact that the corresponding vacuum solutions have no event 
horizons (see eq. (6.22) in Ref. [18J) and are globally nonstatic. 

With no apparent horizons, no black holes may form (more precisely, the whole Universe 
is one black hole). This excludes the quasi-plane and quasi-hyperbolic models from an 
important area of application of the quasi-spherical and spherically symmetric dust models. 



7 The toroidal plane symmetric model in the Szekeres 
coordinates 

The metric of a torus given by constant t and z in (13.21) is ds2 2 = R 2 (dp 2 + dg 2 ), where 
2nR is the circumference of the torus in the p-direction and in the g-direction. In the 
following, we will consider the torus with R = 1, and we will call it the 'elementary square' 
or 'elementary torus'. It will be more convenient to assume that, in the coordinates of 
(13.21) . the elementary torus is the square {p, q} G [— tt, tt] x [— tt, tt], shown in Fig. HI rather 
than {p,q} G [0, 2tt] x [0, 27r]. The IR 2 -space of the (p, q) coordinates can be imagined as 
filled with infinitely many copies of this square. 

As observed in Ref. [18] , the function S in the quasi-plane model can be absorbed into 
the other functions by the redefinition 

(R, E, M) = (R/S, E/S 2 , M/S 3 ), (7.1) 

so we can assume S = 1 with no loss of generality. We do this in the following. 

The coordinates of points to be identified are related in a more complicated way in the 
Szekeres coordinates of (12. ip . in which the plane symmetric model is given by: 

,9 ,9 R,z 2 , 9 „9 4(dx 2 + dw 2 ) . , 

ds 2 = dt 2 ^frdz 2 - R 2 ^ ^ 7.2 

2E ^ [{x - Pf + (y- Q) 2 ] 

with P and Q being arbitrary constants. A line q = qo corresponds, in the (x, y)- 
coordinates, to 

(x - P) 2 + (y - Q - l/q ) 2 = 1/V, (7.3) 

which is in general a circle of radius l/q and the center at (x, y) = (P,Q + l/q ). In the 
special case q = the image becomes the straight line y = Q. Consequently, the lines 
q = ±7r in the (p, g)-coordinates go over into the circles 

( x -Pf + (y-Q T l/ 7 r) 2 = l/n 2 , (7.4) 

while the lines p = ±7r go over into the circles 

( x -P T l/n) 2 + (y-Q) 2 = l/n 2 . (7.5) 
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Figure 4: A map of the elementary torus in the (p, q) coordinates (the central clotted square). Its left 
edge coincides in space with the right edge, the lower edge coincides with the upper edge. Sometimes it is 
convenient to consider the torus as a subset of the M 2 plane, in which case the plane should be imagined 
as covered with an infinite number of copies of the elementary square. It will later turn out that the 
identifications have to be done with a twist, i.e. with reflections in p = and q = 0. The open circles and 
the full circles show the pairs of points to be identified. 
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Figure 5: The image of the torus {p, q} 6 [— tt, tt] x [— n, tt] in the Szekeres coordinates. The values of 
P and Q were chosen arbitrarily, but other elements of the figure are drawn to scale. The small empty 
circles and the small solid circles mark pairs of points to be identified in the general, nonsymmetric case 
(see Sec. EJ. More explanation in the text. 



The image of the central point (p, q) = (0, 0) is the infinity of the (x, y)-p\ane. Conversely, 
the point (x, y) = (P, Q) is the image of the infinity of the (p, g)-coordinates. 

Moreover, from (17. 3ft follows that the image of the area {q 2 < go 2 } (an infinite strip of 
the (p, q) plane contained between q = — g and q = go > 0) is the area outside the circles 
(x — P) 2 + (y — Qt 1/<?o) 2 = l/<?o 2 - Similarly, the image of the area {p 2 < p 2 } is the area 
outside the circles (x — P =F 1/Po) 2 + {y — Q) 2 — 1/Po 2 - Consequently, the image of the 
elementary torus in the (x, y) coordinates will be the infinite subset of the M. 2 plane lying 
outside all four circles, see Fig. 

This is the explanation to Fig. The image of the line p = in the (x, y)-coordinates 
is the vertical line x = P in the figure, with (x, y) = (P, Q) being the image of infinity 
of the (p, q) coordinates. Similarly, the image of q = is the horizontal line y = Q. In 
the (p, g)-coordinates of (13. 2p . the torus is the area encircled by the straight lines p = —tt, 
p = tt, q = —tt and q = tt. The image of the line q = —tt is the circle CI, of radius 1/tt. 
The image of the line q = tt is the circle C2, of the same radius. The image of the torus 
must be contained outside these two circles - in the area covered with vertical strokes. 
Then, the image of p = —tt is the circle C3, and the image of p = tt is the circle C4, both 
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of the same radius l/n. Consequently, the image of the torus must be contained outside 
these circles - in the area covered by horizontal strokes. Thus, the image of the whole 
torus is the common subset of these two areas - the area in the figure outside the thick 
line and covered with crosses. The area inside the thick closed curve contains an infinite 
number of images of copies of the elementary torus. 

Each circle in the figure with the center at (x, y) = (P, Q) and with radius a > 2/n is 
an image of a circle of radius 2/a < it centered at (p, q) = (0, 0) that lies all within a single 
copy of the elementary torus (one such circle is shown in the figure with a dotted line). In 
particular, this applies to a circle of radius 1. 



8 A nonsymmetric toroidal planar model 

We will verify now whether the toroidal interpretation can be applied to the full nonsym- 
metric planar Szekeres model. 

In the plane symmetric subcase, with the metric (13.21) . it would in fact be possible to 
identify points for which p = p Q with points for which p = —p at any fixed po because the 
metric is invariant under the reflections p —>■ —p; and the same is true for q. The discrete 
symmetries in the general planar case, (12.71) with s = 0, are more restricted and they will 
point out more precisely which points can be identified with which. 

In the general planar Szekeres model, as seen from (17.51) . identifying the line p = tt 
with the line p = — tc means identifying the circle C4 in fig. [SJ on which 

{ X -P) a + (y-Q)' = -{x-P), (8.1) 

7T 

with the circle C3, on which 

(x-P) 2 + (y-Q) 2 = --(x-P), (8.2) 

7T 

To see what effect this has on the metric,we write out eq. (12.71) with e = and S — 1: 



ds 2 = dt 2 ' 



2E z) 



nD ( X -P)P, z+ (y-Q)Q, 

,z+ (x-py + (y-Qy 



2 



dz 2 



R 2 (dx 2 + dy 2 ) 

[{x-py + iy-Qyf 

Identifying (18 .ip with (18.21) is clearly possible within the (x, y)-surface, where this is an 
isometry. To make it an isometry also in the spacetime, we have to do this identification 
simultaneously with the substitutions 

x-P^P-x, y-Q^Q-y. (8.4) 

The first is a reflection in the x = P axis in Fig. 0, which clearly must accompany the 
identification of (18.11) with (18. 2p . The second of (18. 4p is a reflection in the y = Q axis, 
which is anyway required if we identify the circles C\ and C^- 

The conclusion is thus: the toroidal interpretation is admissible also in the general case, 
provided the identifications within each (x, y) surface occur in the way marked by the large 
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points in Figs. H] and [5j i.e. with reflections in the x = P and y = Q axes. A square with 
sides identified in this way is in fact a projective plane (for an elementary introduction 
to projective geometry see Ref. [26]). Such a topology of the (x,y) surfaces solves some 
problems, but at the same time creates some others. We will come back to this later. 



9 Shell crossings in the toroidal model 

We showed in Ref. [IB] that the planar Szekeres metric, if its (x, y) surfaces are interpreted 
as infinite planes, always develops shell crossings unless S, z /S is constant, in which case 
the spacetime becomes plane symmetric. The toroidal topology can solve this problem. 

Let us recall, after Ref. [18], that the surfaces of constant t and constant z in the planar 
model, (12. 7p with e = 0, can be parametrised by the variables (9, 0) related to (x, y) by 

2 2 
x-P=-cos0, y-Q = -sm.(j) (9.1) 

(we assumed that S was absorbed into the other functions as in (17. ip ). The value of 1/9 
is the distance from the point (x, y) = (P, Q) in Fig. [51 A shell crossing occurs at such 
(x, y, z) where R, z —R£, z /£ = in (12. 7p . while E ^ 0. In the (9, 0) variables we have 

= -#(P, 2 cos0 + Q, 2 sin0), (9.2) 

and so the locus of a shell crossing is at 

q = thl . (9.3) 

R (P, 2 cos + Q, z sin 0) 

As long as (P, 2 cos + Q, z sin 0) is not identically zero, i.e. as long as P, z and Q, z do not 
vanish simultaneously, this equation has a solution for 9$ However, if the 9 implied by 
(19.31) is large enough, then the corresponding 1/9 will be small, and the point with such a 
value of 9 will lie inside the closed curve in Fig. [5J i.e. outside the elementary torus. For 
this to be the sufficient condition is that the 1/9 of (19. 3p is smaller than y/2/ir - 

the smallest distance of the contour in Fig. [5] from the origin at (x, y) = (P, Q)\ 



R (P, 2 cos + Q, z sin t 



R, z 



< A (9.4) 

7T 



Since |P, z cos0 + Q, 2 sin0| < a/ P, 2 2 + Q, z 2 , a sufficient condition for (19. 4p to hold at all 
values of is 

op 2 

PJ + QJ<^ (9.5) 

at all t and all z. Now we have to verify whether this condition can hold during all of the 
evolution with P, 2 2 + Q, z 2 > 0. We know (see, for example, Refs. [27], [15] or [16]) that 
P, 2 can be achieved at all times and all z - this is identical to one of the conditions 
that guarantee no shell crossings in the Lemaitre - Tolman model. However, we have to 
verify that (19. 5p can still be fulfilled with P, 2 2 + Q, z 2 > when R — > oo. 



3 By definition, 9 > 0. But the correct sign of 9 can always be achieved by changing 
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For this purpose, we follow Hellaby and Lake [27] and use the formulae of Ref. 
Section 18.10. We do not have to consider the case when E = in (12. 8p because, as stated 
earlier, in the planar model E can vanish at isolated values of z, but not on open intervals 
(in fact, E — > only at the asymptotic origin). Thus, the behaviour of the model at E — > 
can be calculated from the formulae given below in the appropriate limit. With E > and 
A = 0, the solution of (12.81) can be represented in a parametric form as 

M 

R(t,z) = ^ (cosh 77 - 1), 
sinhry-r/ = L-j— [ t - t B (z)] , (9.6) 
where t B (z) is the (arbitrary) bang time function. From here we then find 

R M E \2 J M B ' z iyiJ 



(9.7) 



where 

def sinh^sinhrz-r/) ^ « smlir? (9 . 8) 

(coshr/ — l) 2 ' ' (coshr/ — l) 2 

The functions $3(77) and $4(77) have the following properties: 

lim^o $a(v) = 2 / 3 - lim $ s(r/) = 1, 

lim^o *4(r/) = 00, lim $4(77) = 0, 

d$3 d$4 

— - > for 77 > 0, — - < for n > 0. (9.9) 

dr/ dr? 

Then, as shown for the Lemaitre - Tolman model, R, z > for all (t, z) if and only if 

M, 2 >0, E, z >0, t B , z <0. (9.10) 

These can be fulfilled by the appropriate choice of the arbitrary functions. Consequently, 
the properties of the functions $3(77) and $4(77) guarantee that R, z / R > for all t and z, 
and so, with sufficiently small nonzero \P, Z \ and \Q, Z |, eq. (19.51) can always be fulfilled, 
i.e. shell crossings can be avoided throughout the evolution. 

This conclusion remains true when we require that 1/9 < iV2/ tc in (19.31) - (19 .4p . where 
i is a natural number, i.e. that the region free of shell crossings is a multiple of the 
elementary square. This observation will prove useful below. 



10 Formation of structures in the planar model 

As shown for the Lemaitre - Tolman models (see Refs. [25] and [16] . Sec. 18.19), in the 
ever-expanding case E > an increasing density perturbation, p, z / p, freezes asymptoti- 
cally into the background - i.e. it tends to a finite value determined uniquely by the initial 
conditions. Consequently, it is impossible in these models to describe the formation of 
condensations that collapse to a very high density, such as a galaxy with a central black 
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hole. Since the evolution of the quasi-plane and quasi-hyperbolic models is described by 
the same equations, they will suffer from the same problem. (And we have already found 
in Sec. [6] that these models cannot describe black holes.) 

Thus, these models can be used for considering the formation of moderate-amplitude 
condensations and voids. 



11 Interpretation of the mass function M(r) in the 
toroidal quasi-plane model 

Let us recall some basic facts about the mass function M(r) in the quasi-spherical model. 

The amount of rest mass within a sphere of coordinate radius z at coordinate time 
t, when z = z c is the center of the sphere, equals A4 = L p^/\g 3 \d^x, where V is the 
volume of the sphere, p is the mass density given by (12.111) and g% is the determinant of 
the 3- metric t = const in (12.71) . Substituting for p and for g% we get 



4:71 



dx 



dy / du 



M 1U (u) _ 3M£, U 
_VT+2E£ 2 VTT2ES 3 



11.1) 



where u is the running value of z under the integral. Note that £ is the only quantity that 
depends on x and y, it is an explicitly given function, and so the integration over x and y 
can be carried out. As is easy to verify, we have0 



dx 



' 1 
dy^ 2 



47T, 



dx 



dq -£^ 



0. 



'11.2) 



(The first of these just confirms that this is the surface area of a unit sphere.) Using this 
in (lll.ip we get 

r m 

M= / ; (u)du, (11.3) 

VTT2E ' v 1 



which is the same relation as in the LT model, and shows that l/\/T+ 2E is the relativistic 
energy defect /excess function (when IE > and 2E < respectively). 

If we interpret the quasi-plane model as having the toroidal topology, then the calcu- 
lation given above can be done in a very similar way. The integrals over x and y in the 
analogue of (111.11) should now extend over the area S outside the clover-like contour in Fig. 
[5j Just as in the quasi-spherical case, since only £ depends on x and y, we can integrate 
over x and y first. For this purpose, it will be convenient to transform (x, y) by (13.11) . since 
this will greatly simplify the definition of the area of integration: in the (p, q) variables, 
this is just the square [— x [— 7r, 7r]. We have: 



d 2 xy 



(p 2 + q 2 



11.4) 
11.5) 
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This fact was first proven by Szekeres [2]. 
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The volume-element in calculating the mass is non-oriented, so in transforming the integral 
we take the absolute value of the surface-element. Thus: 



d 2 xy— 
d 2 xy-^ 



dp dq = Air 2 , 

J —TT 

dp / dg (pP, z +qQ, z ) = 0. 



11.6) 
[11.7) 



However, here we encounter a problem with the geometrical interpretation. A projective 
plane P 2 is one-sided, and so non-orientable [26] . Thus, it is impossible to divide a 3-space 
in which P 2 is immersed into the interior of P 2 and its exterior. We will clearly need such 
a distinction if we want to speak about a mass contained inside P 2 . (The surface element 
f lX 1 .41) gives no warning about the change of orientation after going all around the plane, 
but this must be an artifact of the coordinate map.) This problem may be solved, for 
example, as follows. In Fig. HI we identify the line p = n with the line p = —tt not in 
the same elementary square Di, but in an identical duplicate copy \3 2 of the elementary 
square. After proceeding from p = —tt to p = tt in n 2 , we identify the set p = tt of \3 2 with 
p = —tt of □].. Then we repeat the trick in the g-direction: we identify the line q = tt in 
□1 with q = —7i in D 3 , still another duplicate copy of 0\. In D 3 , we identify the set q = tt 
with q = —71 of Di, and the set p = tt with p = —71 in D 4 , a fourth copy of □].. The setup 
is illustrated in Fig. [6j In this way, while going across the elementary square, we can arrive 
back to the starting point after crossing at least twice its length. The surface obtained 
in this way is two-sided and orientable. Its total surface area is then 16tt 2 instead of the 
4tt 2 calculated in (jll.fip . but still finite, and the integration in ( 111.61) effectively covers the 
region E four times, so we denote the area of integration by 4E. 

Consequently, the analogue of (jll.ip for e = is 
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As already stated (see the comment around eq. (15.101) ). if the regularity conditions are 
fulfilled, the integral above is finite, and Ai(z c ) = 0. Consequently, we can interpret M(z) 
as the active gravitational mass within a volume extending from the origin at z = z c to 
the running value of z, and 1/V2E is the value of the mass defect /excess factor. 

With the toroidal interpretation, the proof that M(z) is a measure of the active gravi- 
tational mass was rather simple. The same may be shown also when the quasi-plane model 
is interpreted as infinite in extent, but in a more complicated way. We do show it in the 
next section - however, this is only for mathematical completeness. As demonstrated in 
Ref . [IE] , with the infinite spaces the quasi-plane Szekeres solutions have irremovable shell 
crossings, and so are in fact not acceptable as cosmological models. 



12 Interpretation of M(z) for the quasi-plane model 
with infinite spaces 

The consideration below was devised for the model with infinite spaces of constant t, but 
in fact it applies just as well with the toroidal spaces. 
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Figure 6: Four identical copies of a projective plane patched together to form a two-sided surface (the 
4S of eq. (|11.8[) ). Numbers in circles identify the edges of the elementary square; Roman numbers label 
the copies of the elementary square. Copy II is attached to copy I so that it mimics the twist in the 
projective plane. Similarly, copy III is attached to copy I in the same way. Copy IV is attached to copies 
II and III to close the picture. The right edge of the large square is then identified with the left edge 
without a twist, its upper edge is identified with the lower edge also without a twist. The surface thus 
obtained is two-sided, so orientable, and it makes sense to speak about the volume inside it. 



Recall that with e = we are free to rename the functions R, E and M as in ( 17. II) . and 
the result will be the same as if S = 1. Thus, we assume S = 1 throughout this section. 

For the beginning let us consider the plane symmetric subcase of the e = model, which 
has P, z = Q, z = 0. With S = 1, we can write an analogue of the integrals in (lll.2p . Let us 
choose a circle of radius 1 centred at (x, y) = (P, Q) (both P and Q being now constant) 
in every surface of constant t and z. Let d 2 xy be the surface element in the (x, y) plane, 
and let U be the outside of the unit circle. This region has a finite surface area. (From the 
considerations of Sec. [7] we know that if we interpret the (x, y) surfaces as tori then the 
region U will all be contained within the elementary torus.) Then, introducing (u,(p) by 



x 



P + u cos <p, y = Q + usimp, 



we get 
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in every (t = const, z = const) surface. Now let V be a 3-dimensional set in a t = const 
space, extending from z = z c to a running value of z, whose every section of constant z 
is U - the outside of the unit circle [x — P) 2 + (y — Q) 2 = 1. Calculating, as in (111.11) . 
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Thus, in this case, M behaves as the active gravitational mass contained outside a tube of 
coordinate height (z — z c ) which has radius equal to 1 at every z value, while 1 / \JlE plays 
the role of the mass defect/excess factor. We recall that the (x, y) coordinates of (12.71) . 
in the plane symmetric case e = = £, z , S = 1 are related by the inversion (13.11) to the 
Cartesian coordinates (p, q) in a plane, so the outside of the tube in the (x, y) coordinates 
is in reality the inside of the same tube in the Cartesian coordinates. Thus, physically, M 
is the active mass within a tube V 2 - 

We calculated the integrals in H12.2[) and (112.31) around the central point (x, y) = (P, Q). 
However, with plane symmetry, the origin of the Cartesian coordinates can be transferred 
to any other point by a symmetry transformation!^] 

Let us consider the transformation (13. ip . after which the metric becomes (13.21) . which 
is formally (12. 7p with £ = 1 and (x,y) renamed to (p,q). In this form, the transformation 

p = p' + A p , q = q' + A q (12.4) 

(with A p and A q being arbitrary constants) is an isometry of (13.21) . Thus, the transfor- 
mation (112.41) does not change either the metric (13. 2p or the value of the integral (112.21) . 
which, in the variables (p, q), becomes simply 4 J s dpdq = Att, independently of where the 
centre of the circle Si is located. 

Now let us consider the planar metric e = that is not plane symmetric, i.e. with P, z 
and Q, z not vanishing simultaneously. Let U(z) be the outside of a unit circle in an z = 
const surface, with the centre at (x,y) = (P(z),Q(z)). Within each single such surface, 
applying the transformation of variables (I12.ip . we get 

I /*27r /"oo 



d 2 xy — = / dip —du = Att, 
£ Jo Ji u 



£, z f 2n f°° -4u cos <fP, z -4u sin (pQ 



d 2 xy j§ = J o ^ J inr^iziif, — du = (12 . 5) 

These integrals do not depend on P, z or Q, z , but the centres of the circles no longer have 
the same (x, y) coordinates at each z. Thus, to use (112.51) in an analogue of (112. 3p . we 
have to take a volume V which is a wiggly tube: its every cross-section with a constant 
z surface is a unit circle, but the centres of the circles do not lie on a line orthogonal to 
the z = const surfaces. Instead, they lie on the curve in the t = const space given by 
the parametric equations x = P(z),y = Q(z). Because of the second of (112. 5p . (112.31) still 
follows for this single tube. 

The whole 3-space t = const is now no longer homogeneous with respect to the group 
of plane symmetries. However, each single z = const surface in that space is homogeneous. 
In particular, the surface containing the base of the tube, z = Zq, is homogeneous. Thus, 
we can apply the inversion (13.11) with P = P(zq), Q = Q(zq). The inside and outside of the 
unit circle in the z = Zq will thereby simply interchange, but the resulting transformations 
in other z = const surfaces will be more complicated, and the wiggly tube will deform 



n The circle of unit radius in the Cartesian coordinates, when moved to another point of the (p, q) plane, 
will not have a unit coordinate radius in the (x, y) coordinates, and the image of the center of the circle will 
not be the center of the image circle. However, the surface area of the circle and the invariant distances 
between points are not changed. 
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substantially. Still, in the inverted coordinates we are now free to move the centre of the 
base circle (within the z = Zq surface) to any other point. 
We now carry out this plan. Let us denote: 

P(z ) d = f P , Q(^o) = Qo, 

V = (P -P) 2 + (Qo-Q) 2 . (12.6) 

To the metric (12.71) with e = and S — 1 we apply the inversion adapted to the surface 
{t = const, z = z }: 

V _ n q 

z~, 2' 1/ — VO H 5— 

pz _|_ gz pz _|_ gz 

After this, the 2- metric R 2 (dx 2 + dy 2 ) /E 2 becomes: 



x = p o + ^r^i y = Qo + zr L -3- ( 12 - 7 ) 



ds 2 2 = i (dx 2 + dy 2 ) , 

2£ = V(p 2 + q 2 ) + 2(P -P)p + 2(Q -Q)q + 1. (12.8) 

In these coordinates, the surface {t = const, z = z } is explicitly homogeneous, so we are 
now free to shift the origin of coordinates to any other point by 

p = p' + A 1 , q = q ' + A 2 , (12.9) 

with Ai and A 2 being constants. After the shift, the metric is still Szekeres with e = 0, 
but with complicated expressions for the new P, Q and S. 

After the two transformations, (112.7ft and (112.9j) . the region U of integration in (112.51) 
(which was the outside of a tube extending out to infinity) goes over into a finite region - 
the inside of a certain tube whose edge is the image of the family of unit circles in (x,y). 
In the integrals (112. 5p . the two transformations are just changes of integration variables, 
so the values of the integrals do not change, and thus (112.31) still applies. This shows that 
over each point of the surface of constant t and of z = Zq in the Szekeres e = metric, 
we can find a region of finite volume (a wiggly tube) such that the function M can be 
interpreted as the active gravitational mass within that tube. 

As an illustration of the considerations about wiggly columns and their transport over 
the manifold, we now consider a special case of the transformation (112.91) with A\ == A 
and A 2 = 0. But first we give the complete transformation that will take us back to the 
coordinates of (12.71) in the surface z = z . 

To the variables (p f , q') of (112. 9p we apply the inversion in a circle of radius 1 centred 
at (p',q r ) = (0,0), and the shift by (Po,Qo) to the resulting (x',y') coordinates. Calling 
the final coordinates (x 2 ,y 2 ), we calculate the effect of (112.91) on the variables (x,y). The 
complete transformation from (x,y) to {x 2 ,y 2 ) is 

x 2 = P + (x - P + A X U) , y 2 = Q + (y - Q + A 2 U) , (12.10) 



where 



U = (x-P ) 2 + (y-Q c 
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w 



def 



1 + 2 [A x {x - P ) + A 2 (y - Qo)} + ( V + A 2 2 ) U. 



:i2.11) 



This whole set of transformations does not change the metric (12. 7p in the hypersurface 
z = zq, but after the transformations the unit circle U = 1 goes over into the circle 



(x 2 - P - ^i/7) 2 + (V2 -Qo- A 2 hf = l/ 7 2 



;i2.i2) 



where 7 d ^ f A x 2 + A 2 2 - 1. 

Now we specialize this to the 1-parameter subgroup A\ = A, A 2 = 0, i.e. to the shift 
along the p-direction in ( 112.91) . The transformation (112.101) - (112.111) becomes: 



x 2 



W 



p ^ x-P + \{(x-P ) 2 + (y-Q ) 2 } ^y-Qo 

Po H 777 > 2/2 = Qo + 



del' 



w 

\ 2 {(x - P ) 2 + (y- Qo) 2 } + 2X(x - P ) + 1, 



W 



;i2.13) 



and its inverse is obtained by replacing A with (—A), i.e. 

x 2 - Po - A{(x 2 - P ) 2 + (V2 - Qo) 2 } 



x = P + 

def 



Wa 



y = Qo + 



V2 - Qo 

Wo 



Wo = \ 2 {(x 2 - Po) 2 + (y 2 - Q ) 2 } - 2\(x 2 - P ) + 1. 
The Jacobians of the two transformations are, respectively 

J -- 

J -- 



[12.U) 



d(x 2 , y 2 ) 
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d(x,y) 




d(x,y) 
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d(x 2 ,y 2 ) 


= W 2 



;i2.i5) 



The following identities are useful in calculaitons: 



(x 2 - Po) 2 + {y 2 - Qo 
(x-P ) 2 + (y-Qo 



(x - Pq) 2 + (y- Qpf 
Wo 

(x 2 -P ) 2 ^+ {y 2 -Q c 

Wo 



;i2.i6) 



The transformation (112. 13[) takes the circle (x — P) 2 + (y — Q) 2 = u 2 to a shifted circle 
with a different radius, namely 



(x 2 - A) 2 + (y 2 - Bf = (u/po) 2 , 



;i2.i7) 



where 

B 
Po 



def A 2 P [(Pq - Pf + (Qq - Q) 2 -u 2 }-\ \u 2 + P 2 - P 2 - (Qp - Q) 2 ] + P 

Po 

d _e f A 2 Q [(Po - P? + (Qo - Q) 2 - u 2 ] - 2AQ (Po -P) + Q 



def 



Po 

A 2 [(P - P) 2 + (Qo - Q) 2 - u 2 ] - 2A(P - P) + 1 
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'12.18) 



(in fact, these will be applied with u = 1.) 
Now using PH| - (TT2~T6D we find: 



- Pf + (y -Qf = S^ [{x 2 + af + (y 2 + (3) 2 ] , (12.19) 
W 



where 



S ^ 1 - 2A(P - P) + A 2 [(P - Pf + (Qo - Qf] 

P -P-\[(P -py + (Q -Qf] 



def n 
a = -Pn + 



So 



(3 = -g + ^V^- (12.20) 
bo 

Note that and (3 are functions only of z, they do not depend on x 2 and ?/2- We can 

see that the Wq that will appear in the transformed integral in (112. 5p . d 2 x?//£ 2 , will be 

canceled by the Wq from the Jacobian of the transformation, and what remains will be 
1 4 f 1 

— rdxd?/ = ;r / Trdxodyo, (12.21) 

£ 2 So{z) 2 J [( X2 + a f + (y 2 + m 2 

an integral of exactly the form (112. 5ft . except for the additional factor 1/ So 2 ■ 

The transformation (112.131) affects also the metric. Under (112.131) . £ changes as follows: 

£ = ( x - Pf + ( y - Qf = S [{x 2 + af + (y 2 + Pf] . (12.22) 

Wo 

The expression dx 2 + dy 2 , after the transformation ( 112. 131) goes over into 

(dx 2 2 + dyl) /W 2 , (12.23) 
and so the two equations above imply that after the transformation: 



dx 2 + dy 2 dx 2 2 + dy 2 2 



So 2 [(x 2 + ay + (y 2 + pyy 



;i2.24) 



We are in the same Szekeres model as at the beginning, but at a different location. The 
alien form of the Szekeres metric in these coordinates results from the fact that the transfor- 
mation (112. 1 3j) is not a symmetry, apart from the single surface z = zq. Since the quantities 
in this equation (a, f3 and So) depend on the continuous parameter A, the equation above 
describes the result of a shift to any location (with A = corresponding to identity). 
Equations ( 112.241) and ( 112.211) show that the transformation ( 112.131) is area-preserving. 

The second integral in (112.51) does not change its zero value under the transformation 
( 112.131) - since that transformation is simply a change of variables under a definite integral, 
applied both to the integrand and to the area of integration. Thus, the second integral in 
( I12.5P does not contribute to calculating M. in ( 112.31) . 
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13 Summary 



Continuing the research begun in Ref. [18] . geometrical properties of the quasi-plane 
Szekeres model were investigated here, along with the corresponding properties of the 
plane symmetric model. The following results were achieved: 

1. The pattern of decelerated expansion in the plane symmetric model was analysed 
and shown to be in complete disagreement with the Newtonian analogues and intuitions 
(Sections [3] and H]). An example of a Newtonian potential that gives a similar pattern of ex- 
pansion has co-axial parallel cones as its equipotential surfaces; it has not been investigated 
whether such a potential can be generated by any realistic matter distribution. 

2. Embeddings of the constant t and constant z surfaces in the Euclidean space suggest 
that the flat surfaces contained in the plane symmetric model can be interpreted as flat 
tori whose circumferences are proportional to the function R(t, z), and thus vary with time 
(Sec. E]). Such a topology immediately explains the pattern of expansion and implies that 
the total mass contained within a z = const surface is finite. 

3. The quasi-plane and quasi- hyperbolic models are permanently trapped (Sec. [6]), so 
no apparent horizons exist in them. Consequently, these models cannot be used to describe 
the formation of black holes (the whole Universe is one black hole all the time). 

4. The toroidal interpretation can be applied to the full quasi-plane model (Sec. [S]). 
However, the flat surfaces defined by the geometry must then be interpreted as projective 
planes rather than tori. 

5. The toroidal quasi-plane model can be free of shell crossings (Sec. [9]), unlike the 
model with infinite spaces. 

6. For both interpretations (toroidal and infinite), the quasi-plane model cannot de- 
scribe the formation of structures that collapse to very high densities (Sec. \TU\i . since the 
density perturbations tend to finite values in the asymptotic future. 

7. To calculate the mass within a volume, the geometrical interpretation of the flat 
surfaces has to be modified yet again, because the projective plane is one-sided and so 
has no interior. Thus, instead of identifying the opposite edges of the elementary square, 
we introduce 3 extra identical duplicate copies of the elementary square and patch them 
together into a two-sided surface. With the toroidal interpretation thus modified, the mass 
function M(z) is proportional to the active gravitational mass contained within a 'radial' 
coordinate z (Sec. [TTj) . With infinite spaces, this function is proportional to the active 
gravitational mass within a 'wiggly tube' of finite radius (Sec. [T21) . 

Whether the toroidal interpretation is a necessity is still unknown. However, this paper 
demonstrated that with the toroidal topology the quasi-plane model becomes in several 
respects simpler, and, however paradoxical this may sound, more realistic. 

The quasi-plane model with toroidal spaces may be a testing ground for the idea of 
a 'small Universe', proposed by Ellis [28]. A small Universe is one with compact spatial 
sections, in which thus a present observer has already seen several times around the space. 
Several papers were devoted to checking this idea against the observational data (see, for 
example, Refs. [29] - |34j ; a conclusive proof of any nontrivial topology is, unfortunately, 
still lacking). However, the background geometry has always been a homogeneous isotropic 
Robertson - Walker metric with identifications in the underlying manifold. The quasi-plane 
toroidal Szekeres model has a less general topology (identifications in it occur only in two- 
dimensional surfaces, in the ^-direction the space is infinite), but is inhomogeneous and 
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has no symmetry, so might be useful for considering light propagation and comparing the 
mass distribution in the model with the observed images. 
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